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ABSTRACT 
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This paper uses data obtained from the galaxy luminosity function (LF) to 
calculate two types of radial number densities statistics of the galaxy distribution 
^ ; as discussed in Ribeiro (2005), namely the differential density 7 and the integral 

differential density 7*. By applying the theory advanced by Ribeiro and Stoeger 
(2003), which connects the relativistic cosmology number counts with the astro- 
nomically derived LF, the differential number counts dN/dz are extracted from 
the LF and used to calculate both 7 and 7* with various cosmological distance 
Q-f definitions, namely the area distance, luminosity distance, galaxy area distance 

and redshift distance. LF data are taken from the CNOC2 galaxy redshift survey 
' and 7 and 7* are calculated for two cosmological models: Einstein-de Sitter and 

an Q mo = 0.3, ^a — 0.7 standard cosmology. The results confirm the strong 
dependency of both statistics on the distance definition, as predicted in Ribeiro 
^ 1 (2005), as well as showing that plots of 7 and 7* against the luminosity and red- 

shift distances indicate that the CNOC2 galaxy distribution follows a power law 
pattern for redshifts higher than 0.1. These findings bring support to Ribeiro's 
(2005) theoretical proposition that using different cosmological distance mea- 
sures in statistical analyses of galaxy surveys can lead to significant ambiguity in 
drawing conclusions about the behavior of the observed large scale distribution 
of galaxies. 
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1. Introduction 

The galaxy luminosity function (LF) is an important analytical tool for testing cosmo- 
logical models. It conveys observational information on how the local density of galaxies 
varies for different types of galaxies and with differing cosmic environments, as well as how 
the local galaxy density changes with cosmic epoch. The LF has been extensively and sys- 
tematically calculated for a great variety of galaxy redshift surveys, in terms of the number 
of galaxies, their morphological types, observational bandwidths and redshift depths. In 
addition, since in calculating the LF one often includes a model of source evolution, the 
presently available LF parameters obtained from different galaxy samples provide an im- 
portant collection of data for testing cosmological models. To do so, however, we need a 
detailed theoretical framework connecting the astronomical LF data calculation and practice 
to cosmological modeling. This is essential for extracting information that can be used for 
further constraining cosmological models. 

In a recent paper, Ribeiro and Stoeger (2003; from now on RS03) have presented a 
detailed relativistic treatment directly linking the cosmology source counting theory with LF 
data. RS03 started with a general relativistic treatment of number source counts in a general 
spacetime and then derived expressions for relativistic density per source and total relativistic 
energy densities of the universe in terms of the luminosity and selection functions. These 
theoretical/observational relationships were applied to test the LF parameters determined 
from the CNOC2 galaxy redshift survey (Lin et al. 1999) for consistency with the assumed 
Einstein-de Sitter (EdS) cosmological model. Other possible applications and extensions of 
this general framework were also outlined. 

This present work explores one of these extensions of RS03. First, we shall show how 
to extract the differential number counts dN/ dz from LF data. Inasmuch as calculating the 
LF parameters from a galaxy sample requires the assumption of a cosmological model, since 
a volume definition, usually comoving, must be employed, one can factor the cosmology out 
of the LF parameters and re-obtain the "raw" number count data. The methodology for 
determining LF's requires fitting the original data to some analytical form of the LF, usually 
of the Schechter type, where a model of source evolution is also included. Therefore, by 
"'raw' number count data" we mean that the reverse process eliminates the cosmology, but 
does not remove the fitting itself, that is, the LF analytical assumption and source evolution 
corrections. 

Once these "raw" dN/ dz are in our possession, we carry out a statistical analysis with 
the differential density parameters discussed in Ribeiro (2005), which are basically different 
forms of number densities. As such this analysis is based on the differential number counts. 
Thus, it can be done in any cosmological model (since the dN/dz data does not have any 
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kind of cosmology built into it). Here we test the theoretical proposition advanced in Ribeiro 
(2005) that using different cosmological distance measures in statistical analyses of galaxy 
surveys can lead to significant ambiguity in drawing conclusions about the behavior of the 
observed large scale distribution of galaxies. 

Besides the sensitivity of these large scale analyses to the cosmological distances used, 
comparison of number counts at different redshifts is made more difficult by the fact that, as 
we look down our past light cone, we are counting galaxies at different redshifts, and therefore 
at different times. Our line of sight cuts across many different space-like, possibly nearly 
spatially homogeneous, surfaces. This means, that, even aside from the distance measure 
employed, the number counts down our light cone will not correspond in any simple way 
with the integrated galaxy number density over a volume on any space-like three-surface. 

Thus, even though there is "almost homogeneity" 4 of the distribution of galaxies, and 
therefore of mass-energy, on spatial 3-surfaces (Peebles 1993 and references therein, Hogg 
et al. 2005), we should not expect to find "observational homogeneity" as we count galax- 
ies down our past light cone. Yet, employing some distance measures (e.g. the co-moving 
distance and the galaxy area distance) leads to the appearance of such observational homo- 
geneity. For other measures of distance, we obtain, as Ribeiro (2005) has shown theoretically 
and as we here shall show from the CNOC2 data, a power-law galaxy distribution down the 
light cone, indicating a possible self-similar, hierarchical clustering or fractal-like pattern - 
even though the actual spatial distribution is almost homogeneous. This is simply the re- 
sulting effect of collecting data from successively earlier epochs as we go to larger redshifts. 

Besides showing how to extract cosmology-free dN/dz information from LF's, a pri- 
mary objective of the paper is to clarify these issues, emphasizing the illusory and ambigu- 
ous character of such observational homogeneity and of such purely observational galaxy 
distributions, in general. 

The main conclusions we reach in this paper are as follows. The galaxy distribution given 
in the CNOC2 redshift survey supports Ribeiro's (2005) theoretical finding that the statistics 
of the observed galaxy distribution are strongly dependent on the choice of cosmological 
distance definitions. Both the differential density 7, which gives the rate of increase in 
galaxy number counts with distance down the past light cone, and the integral differential 
density 7*, which is simply the integral of 7 over a specified volume, decay with higher 



4 We use this term to describe succinctly one of the key characteristics of a perturbed Friedmann- 
Lemaitre-Robertson- Walker (FLRW) universe. Our universe is obviously not exactly spatially homogeneous 
- otherwise we would not be here! But there is strong evidence that it is "almost homogeneous," or "almost 
FLRW," on very large scales. 
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luminosity and redshift distances, with no flattening of this decrease being found within 
the limits of the CNOC2 survey (0 < z < 1). These decreasing values of the differential 
densities show a tendency to linearization in a log-log plot for z > 0.1, that is, a power law 
pattern, which suggests a self-similar behavior for the observed galaxy distribution in those 
ranges. Again, it is important to recognize that these parameters are measuring changes 
in galaxy number density down our past light cone, not in the space-like surfaces of the 
Friedmann-Lemaitre-Robertson- Walker (FLRW) models where there is evidence of "almost 
homogeneity". These results are valid for both EdS cosmology and the currently favored 
FLRW standard model with £l mo = 0.3 and J1a = 0.7. However, by comparing the results of 
both models it seems that the linearization tendency mentioned above is more pronounced 
in the second model, suggesting that FLRW cosmologies with a less than critical matter 
density parameter might be a better fit as far as a self-similar distribution is concerned. 

The plan of the paper is as follows. In §2 the main equations, derived in RS03, necessary 
for our purposes here are summarized and the method for extracting the differential number 
count dN/dz from LF data is developed and applied to the CNOC2 galaxy redshift survey 
data. In §3 we calculate the differential density statistics in two cosmological models, namely 
EdS and the presently favored standard model with Q mo = 0.3, Qa = 0.7. §4 summarizes 
our conclusions and the appendix briefly describes the error analysis performed with the 
CNOC2 data. 



2. Differential Number Counts from LF Data 

In this section we present a summary of the relevant results from RS03 and derive the 
general equation for the observed differential number count. Then we will specialize these 
equations to the EdS model in order to be able to obtain the "raw" dN/dz from the LF 
parameters calculated with the CNOC2 galaxy redshift survey. 

2.1. General Relations 

We start with the general equation for the number of cosmological sources N in a volume 
section at a certain point down the null cone whose affine parameter distance is given by y, 

dN = (d A f (1 + z) dn n(y)dy. (1) 

Here d A is the area distance (also known as angular diameter distance, observer area distance 
or corrected luminosity distance) of this section from the observer's viewpoint, dflo is the 
solid angle of the volume section at the observer's position, n is the number density of 
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radiating sources per unit proper volume and z is the redshift. This equation is simply a 
more convenient rewriting of Ellis' (1971) key result (see also Plebahski and Krasihski 2006) 
for cosmological number counts in any cosmological model (see details in RS03, §2.1.1). 
The number density n appearing in the equation above is related to the galaxy luminosity 
function <fi as follows, 



n(y) = ^[z 



(p(l)dl, 



(2) 



where / = L/L*, L is the absolute luminosity of the source and is the luminosity scale 
parameter (RS03, eqs. 29, 35). In addition we have that 



p = M„^, 



(3) 



where p is the local density as given by the right-hand side of Einstein's field equations and 
M, g is the average galaxy rest mass (RS03, §2.1.2, eq. 40). Therefore, equation (1) can be 
written as, 



^L = (d A ) 2 (i + z) dn 

dz 



P 



M, 



dy 

dz 



(4) 



This equation can be directly written in terms of observable quantities in a certain 
bandwidth W as follows (RS03, §4.4, eqs. 46, 84), 



(5) 



~dN~ 




— ^ ^ a w 

W 


~dN~ 


w 


dz 





dz 





Here the subscript zero means an observed quantity, and the superscript index w means the 
observational bandwidth used when the data was collected. The factor a w is the fraction of 
galaxies in a certain wave band that are not counted in other bandwidths. This is necessary 
to avoid overcounting due to possible overlap of galaxy counts in different bandwidths. 
Therefore, we have that (RS03, eqs. 45, 46), 



a w (z) 



1, for W = l, 

b w (z) < 1, for W>1. 



(6) 



where, b w (z) is the fraction of galaxies in each waveband, that is, when W > 1, which are 
not counted in wavebands 1,2, ... , (W — 1). Thus, equation (4) can also be written 



dN 

dz 



(d A ) 2 (i + z) dn 



dy 

dz 



P 



(d A ) 2 (i + z)dn ^Yl 



dz 



P 



Similarly, equation (3) yields 



(7) 



(8) 
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Here i\) w is the selection junction calculated from the luminosity function obtained in a 
certain wave band W, 

^[l(z)\ = / 4> W {1) dl, (9) 

is the galactic morphology population fraction, giving the abundance of each galactic 
type relative to the total number of counted galaxies in each redshift range and observed 
in a certain wave band W. Finally, M v is the average rest-mass for galaxies of a certain 
morphological class v (see eqs. 13, 14, 87 of RS03). 

From the discussion above it is clear that it is still possible to write the observed ratio 
between the local density and the average galaxy rest-mass as follows, 



P 



P?M V *I>. 



E V P V W M V 



(10) 



This equation allows us to write the observed differential number counts in terms of the 
chosen cosmological model and observational quantities extractable from galaxy catalogues. 
This result may be written as below, 



dN 

dz 



{d A f 



'1 + z) dQ.Q — > a, 



dy 

dz^' 



E v K W M V 



;n) 



Notice again that this equation is absolutely general, independent of cosmological model and 
independent of the survey, since its observational part under the summation signs is written 
in such a way as to be applicable, in principle, to any galaxy sample (RS03, §3.1). 



2.2. Evolution of Galaxy Masses 

A straightforward examination of equation (11) shows that it depends on the typical 
average rest-mass Ai v of galaxies belonging to the morphological class v. The way A4 V 
appears in this equation might indicate that this typical mass remains unchanged for different 
redshifts, a hypothesis which is certainly unrealistic as there is today evidence for galaxy 
mergers. Therefore, we must consider some dependency of the typical galaxy mass with the 
redshift in order to have a realistic assumption. 

One way of doing this would be by replacing M. v in equation (11) by [f(z)A4 v ], where 
M. v is the typical mass of v type galaxy at z — 0, that is, at rest with us, or here and now, 
and such that the function f(z) obeys /(0) = 1. So, if a galaxy doubles its mass from z — 1 
to z = we could take f(z) = (1 + z)' 1 . A milder growth would be achieved by assuming 
f(z) = (1 + z)~ q , with < q < 1. So, for q = 0.5 the galaxy would have 70% of its present 
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mass when z — 1, whereas with q = 0.2 we go from 0.87 M v at z — 1 to the present (z = 0) 
rest-mass value M. v . 

There is support for this kind of galaxy mass evolution as Carlberg et al. (2000) showed 
that in the CNOC2 and in the Caltech Faint Galaxy Survey galaxies of mass within an order 
of magnitude of AA* will grow by about 0.15.M* from z — 1 to z — 0. Here M* is the typical 
mass of a galaxy with luminosity given by the luminosity scale parameter of Schechter's 
LF (eq. 2). Simulations carried out by Murali et al. (2001) agree with Carlberg et al. (2000). 
Both studies emphasize, however, that these results only pertain to galaxies which are large, 
that is, with masses near Ai*, leaving out low mass galaxies which should also grow by 
mergers and accretion, but not as much as the larger galaxies. Both studies do not mention 
possible dependency of galaxy mass evolution with morphological types. Therefore, a mild 
growth such as setting q = 0.2 for all galaxy morphologies seems a quite reasonable first 
order approximation considering all these uncertainties. 

At this point one must look again at equation (11) to notice that if we replace Ai v by 
[f(z)M v ] the mass evolution estimate f(z) is canceled out in that equation. This means that 
we can use the values of the typical masses at rest with us in equation (11), since, in that 
expression, the following approximation holds: M v ~ M v . In other words, the observed 
differential number counts (11) constructed from LF parameters, which is the basic quantity 
to be used in all analyses of this paper, is, therefore, not affected by galaxy mass evolution, 
at least to first-order approximation. 

Notice that although galaxy mass evolution does not directly change equation (11), 
it appears indirectly since the selection functions in equation (11) were obtained from LF 
parameters whose derivation usually include some sort of source evolution. 



2.3. Einstein-de Sitter Cosmology 

We now specialize the equations above to Einstein-de Sitter (EdS) cosmology. This is 
done by finding the equations for the affine parameter and the area distance in terms of the 
redshift. It is simple to show that both equations can be written (RS03, §4.2) 

dy c 



dz H {l + z) 1/2 ' 



(12) 



d A = p 1 + Z ~^ t (13) 
H (i + 2 ) 2 ' 1 } 

where H is the Hubble constant and c is the light speed. Throughout this paper we take 
H = 100 km s" 1 Mpc -1 as this is the value adopted by Lin et al. (1999). We aim here to 
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factor the cosmology out of LF data and since the value of the Hubble constant is embedded 
in LF calculations, we must adopt the same value originally set for the LF calculation. In 
the case of this paper it then must be the same as used by Lin et al. (1999). 

An important issue in LF calculations is the volume definition. Currently observers 
adopt the comoving volume, whereas all equations above assume the proper volume. So, 
we need a volume conversion, obviously dependent on the cosmological model, in order to 
actually use the selection function in our equations above. Thus, in the EdS model the 
selection functions in both these volumes are related by (RS03, footnote 6 and §4.3), 



(14) 



where the left superscripts PR and C denote that the selection function is evaluated with 
the proper and comoving volumes, respectively. 

Bearing in mind the results above and spherical symmetry, we are now in position to 
write the observed differential number counts (11) in terms of the selection function in the 
EdS cosmological model as follows, 



~dN~ 




/36ttc 3 \ 


dz 





V H ) 



7/2 



(15) 



2.4. dN/dz from LF Data of CNOC2 Galaxy Survey 

The CNOC2 sample was obtained in three wave band filters: R c , U, B AB (hereafter, we 
shall for notational ease, abbreviate the latter as "B"). This gives us the possible values for 
the wave band index W = R c , U, B. So, according to equation (6) we have a w = b w and 
since all galaxies appear in all bandwidths, we have that 

b Rc = b B = b u = 1/3. 

This catalogue has morphological indexes v — 1, 2, 3 denoting the three morphological types 
in which the sample population was divided, early, intermediate and late spectral types, 
respectively (Lin et al. 1999). These galaxy morphological types have very approximate 
dynamical average rest masses equal to Ai± = 0.5 x lO n 7W , M.2 = 0.3 x lO 11 ^© and 
M 3 = 0.1 x lO n Al (Sparke and Gallagher 2000, pp. 204 and 264). As discussed in RS03, 
the CNOC2 sample is such that the proportion P v of each galaxy type population is as 
follows (RS03, §3.4, tables 1-3), 

P*° = = p» = 0.29, 
p«c = pv = ps = o.24, 
P^ = pv = pB = o.47. 
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With these results equation (15) can be applied to the CNOC2 sample, yielding, 



dN 



^36ttc 3 ^ ( 



1 + z- VI + z) 



3 



V 3 P W M ( c vb w ) 

Z^,v=l 1 V ■ /V '» lCNOC2y» J 



dz 




(16) 



J o 



(Rc,U,b) 



where the correction 



( CN oc C 2 C) = ™ 2 (°0 



is necessary in order to have correct units. Considering now the CNOC2 figures above, this 
equation can still be changed to, 



where the C and CNOC2 indexes of the selection function in equation (16) were dropped. 
The observed differential number count of the CNOC2 redshift survey can then be finally 
calculated. Tables 1, 2 and 3 below reproduce the numbers calculated and presented in RS03 
with errors calculated here (see §A), which are enough to evaluate equation (17). The results 
are shown in table 4. 

Notice that the value of the Hubble constant in equation (17) must be the same as orig- 
inally used by Lin et al. (1999) in deriving the CNOC2 LF, that is, H Q = 100 km s" 1 Mpc -1 . 

3. Differential Density Statistics with CNOC2 Derived Data 

We shall now present the differential density radial statistics required to analyze the 
galaxy number density distribution and its conceptual implications. Some of this was previ- 
ously discussed in Ribeiro (2001, 2005) and, therefore, §3.1 is in part a review. Subsequently 
we will apply these general expressions to the CNOC2 data obtained in the previous section 
in two different standard cosmological models. 




x (^f c + f{ + tpf ) 1.45 x 10 10 + 
+ (V^ c + V£ + V> 2 S ) 7.2 x 10 9 + 



(17) 
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3.1. Basic Equations and Concepts 

3.1.1. Definitions 

The differential density 7 at a certain observed distance d is defined by the following 
expression, 

1 dN 

7= W (18) 

where 

So = Mdof (19) 

is the area of the observed shell of radius do- This definition appears to have been initially 
used in the context of the galaxy distribution by Wertz (1970, 1971). As can be easily seen, 
7 is really very much like the similar differential number counts parameter dN/dz (see eq. 
4), except that it is differential with respect to an observable distance down our past light 
cone, instead of with respect to z (which is also an observable parameter down our past 
light cone). Furthermore, it is "normalized" with respect to the "observed surface area" 
So- Essentially, then, 7 tells us the rate of growth in the number counts, or more exactly 
in their density, as one moves down the past light cone in the observable distance do- As 
is obvious, the behavior of this parameter will depend heavily on which distance we employ 
for d - observer area distance, galaxy area distance, luminosity distance, redshift distance, 
or something else. 

Furthermore, Ribeiro (2005) has advanced an extension of 7 called the integral differ- 
ential density 7*, which is simply the integration of 7 over the observed volume V . Such a 
definition yields, 

7 * = _L f 7 dV , (20) 

V JVo 

where 

Vo = gvr(do) 3 . (21) 

We should carefully note that dVo is different from the volume element given in the usual 
form of the metric. This is because it involves d(d ) down our past light cone, instead of 
something like a(t)dr at one given time t — t± - that is on space-like surface, usually one of 
"almost homogeneity." 
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lt is more useful to write 7 and 7* in terms of the redshift, as follows, 5 



7*(*) 



ciiV 
dz 



S {z) ^j- [do(z)\ 



dV 



, (z) ^dz. 

Vn(z) .L dz 



(22) 
(23) 



The equations above can be obtained from any given cosmological model, since all quantities 
on the right hand side can be calculated from the geometry of the model and the assumed 
matter distribution. The equivalent equations built from data derived galaxy samples should 
then read as follows, 



t5(*) 



~dN~ 






-1 


dz 





So ^o)_ 





V Jo 70 



dVo 

dz 



dz. 



(24) 
(25) 



Notice that by their definitions in equations (18), (19), (20), (21), the following result holds, 

7* = ^ = <n>, (26) 
where (n) is the average number density. 



3.1.2. Conceptual Implications 

Spatial and Observational Homogeneities 

In order to interpret the meaning and significance of 7 and 7* further, we emphasize that in 
cosmological models we can define two distinct types of homogeneity. Spatial homogeneity 
is defined in the usual sense, that is, the usual spatial density p appearing on the right 
hand side of Einstein's field equations is dependent only on time, p = p(t), which means 
that in each constant time hypersurface the local density remains unchanged. So, spatial 
homogeneity means that quantities such as density are constant on a space-like foliation of 
surfaces. The primary example, of course, are the surfaces of homogeneity in FLRW. 



5 This notation for the differential densities is slightly different than in Ribeiro (2005), since here we follow 
RS03 very closely and adopt the subscript index zero to refer to observed quantities, built from galaxy sample 
derived observations as much as possible, whereas the unindexed quantities refer to theoretical expressions, 
built from theory. 
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However, there is another type of homogeneity, observational homogeneity, that can 
also be defined in cosmology. It is the property that some average density calculated with 
observational data derived from galaxy redshift surveys is constant. This is sometimes what 
one has in mind when one talks about the possible homogeneity of the galaxy distribution 
in the context of observational cosmology. So, observational homogeneity is defined in terms 
of average observational quantities along our past like cone, instead of in a given space-like 
surface. Operationally, the observational average density (p ) is constant. 6 How to relate 
(p ) to the above-mentioned pit), how to determine pit) from an observed (po), and what 
the relationship between observational homogeneity and spatial homogeneity are key issues. 

Since a basic method for observationally dealing with galaxies is by counting them, it 
is reasonable to state formally that observational homogeneity is achieved when the average 
number density (no) remains constant, as one averages over different volumes down our past 
light cone. So, 

. . N . , 

(n ) = — = constant, (27) 

where Nq is the number of observed galaxies counted within a certain volume Vq. As indicated 
above, this definition of observational density implies that both N and Vo are obtained along 
our past light cone. So, equation (27) is not applicable over space-like surfaces. 

The practical calculation of (no) by means of equation (27) encounters some basic oper- 
ational difficulties. We immediately face the problem of identifying the volume Vo- Actually, 
since N is observed, V should also be an observed quantity. But, we cannot yet translate 
the observed redshift z into a cosmological distance without a definite cosmological model. 
So, we must choose a certain distance measure among several possibilities from the model 
(many workers often choose the comoving distance), observe N and calculate (n ). There- 
fore, although spatial homogeneity can be uniquely defined relative to a space-like slicing, 
observational homogeneity cannot be. Obviously, a model can satisfy spatial homogeneity, 
without satisfying observational homogeneity, or vice-versa. In fact this will often be the 
case. 

This was shown to be correct in Ribeiro (2001, 2005), with earlier discussions in Ribeiro 
(1992, 1993, 1995). Since we can also obtain iV from theory, Ribeiro (2005) showed that, 
although EdS cosmology is spatially homogeneous by definition, it may or may not be ob- 
servationally homogeneous, depending on the particular type of distance chosen. Thus, 
according to Ribeiro (2005), an FLRW model can satisfy both types of homogeneity, spatial 



6 The notation may be a bit confusing here, as po is the local density at the present constant time to, 
whereas (po) is the average of p[t(r)] along the null cone t(r), since /o[t(r)] will change along the past null 
cone (see Ribeiro 2001 for in depth discussions of this topic). 
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and observational, if, and only if one adopts the galaxy area distance, or comoving distance, 
in the definition of observational homogeneity, as given by equation (27). Furthermore, it is 
clear that lack of observational homogeneity is perfectly consistent with spatial homogeneity, 
and that there will be cases of observational homogeneity in which there is no spatial ho- 
mogeneity (Rangel Lemos and Ribeiro 2006, in preparation). In specifying the link between 
the two the distance being used is crucial! 

From this point of view, interpreting 7 and 7* is straightforward. Considering equations 
(26) and (27) it is clear that 7q = (n ). This means that 7 is a useful tool for determining 
(no). But this is true only if we are trying to determine whether or not the galaxy number 
count N is observationally homogeneous. 

Other Remarks 

It is important to point out that only one of the two factors on the right hand side of 
equation (24), can, at present, be calculated directly from observations. Only the observed 
differential count [dN/dz] can be evaluated from equation (11) or one of its specialized 
forms (15), (16), (17). 7 To calculate the second factor of equation (24), the geometrical part 
involving So and d , one must still assume a cosmological model, inasmuch as obtaining 
observational distances, either by standard candles or standard rods, of every galaxy in a 
redshift survey is currently beyond our observational capabilities. 

Along with this, as we have already indicated above and in the introduction, there is 
the key problem that the observational distance do as obtained from theory is not unique, 
meaning that the observed density of the galaxy distribution is dependent on this distance 
choice. This applies even both to 7 and to 7* - 7* depends on the volume Vq over which the 
galaxies are counted, which in turn depends on the distance choice. We have already seen 
that that this choice can make a significant difference. Which distance is appropriate? 

As it has been extensively argued elsewhere, these are the basic, but, so far, not gen- 
erally acknowledged, dilemmas underlying the observational determination of the possible 
smoothness of the universe (Ribeiro 2001). Bearing these points in mind, it becomes clear 
that comparing 7 with 70 and 7* with 7q does not provide a true check of the cosmological 
model against observations, since the assumed cosmological model appears in both pairs of 



7 The observed differential number count [dN/dz] does not need to be estimated only from LF data, as 
it can be directly observed. Nevertheless, in the context of this paper, we are recovering it from the LF's, 
since those are often what is made available by the observers. 
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equations. 8 We cannot yet determine the geometrical factors purely observationally - they 
must be calculated from the cosmological model which we adopt on other evidential grounds. 
Therefore, such a comparison is merely a check of consistency, the best one can do at the 
moment considering the present limitations in determining the characteristics and distances 
of galaxies. 



3.1.3. Cosmological Distances and the Reciprocity Law 

As stated above the observed distance d$ is not unique. There are at least four distinct 
possibilities (Ellis 1971; Ribeiro 2001, 2002, 2005): the observer area distance d A (which is 
equal to the angular diameter distance), already introduced in §2.1 above, the luminosity 
distance d L , the galaxy area distance d G (also known as effective distance, angular size dis- 
tance, transverse comoving distance or proper motion distance) and the redshift distance d z . 
These distances have expressions dependent on the assumed cosmological model. In FLRW 
cosmologies the last one may be written as, 

d, = (28) 

The other three distances are linked to each other by the remarkable reciprocity theorem or 
Etherington's reciprocity law, valid for any cosmological model and which reads as follows 
(Etherington 1933; Ellis 1971; Plebahski and Krasihski 2006, p. 256), 

d L = (l + zfd A = (l + z) d G . (29) 

These distances can, in principle, be observationally determined, since this law requires 
that source and observer should only be connected by null geodesies (Ellis 1971). Fur- 
thermore, if we assume a cosmological model, theoretically defined distances such as the 
comoving distance, the proper distance, the distance defined by Mattig's formula, etc, can 
be reduced to one of the four above (see, e.g., Ribeiro 2001, 2005). 



8 This is also true for any density constructed with galaxy samples (Ribeiro 2001; Abdalla, Mohayaee 
and Ribeiro 2001). 
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3.2. Einstein-de Sitter Model 

In EdS, besides equation (28) the other observational distances are straightforwardly 
calculated in terms of z as being given by (Ribeiro 2001, 2005), 



A 2C 



-(l + z- VTTz) , (30) 

V ' 



and 9 these results allow us to derive expressions for the observed area and volume for each 
distance, i.e., S L , S A , S G , S z and V L , V A , V G , V z in terms of the redshift (Ribeiro 2005). 
Thus we are able to use the results of table 4 to calculate the observed differential density 
(24) and, after numerical interpolations and quadratures, the observed integral differential 
density (25). 

Table 5 shows the results for both differential densities in each observational distance, as 
well as the values for each distance in the redshift range of the CNOC2 sample. The figures 
show clearly that the general conclusions about the dependence on the distance choice in 
measuring the density of the distribution of galaxies hold. Both the observed differential 
density 70 and the integral differential density 7o increase for higher z when written in terms 
of the area distance d A , whereas both decrease when they are written in terms of d L or d z . 
There is a slight increase when d G is used, but that is mostly within the error margins and, 
even so, so small for both [7 G ] and [7*] that one can consider that in practice both density 
statistics remain unchanged when written as functions of the galaxy area distance d G . Since 
this last distance is the same as the comoving distance (Ribeiro 2005), widely used in testing 
both the possible homogeneity or inhomogeneity of the galaxy distribution, it seems that the 
conclusion reached in Ribeiro (2001) and (2005) about d G being inappropriate for measuring 
the possible inhomogeneity of the Universe do seem to hold, at least as far as the CNOC2 
survey is concerned. Comoving distance is only suitable to test spatial homogeneity in FLRW 
models. These results are shown graphically in figures 1 and 2. 

Figures 3 and 4 show both differential densities obtained from the CNOC2 data and 
plotted against the luminosity and redshift distances, respectively. 7 and 7* decrease as z in- 
creases, showing a tendency towards a power law behavior at scales of z > 0.1. Since power 



In EdS model the comoving distance is equal to d G , apart from a constant (Ribeiro 2005). 
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laws are usually related to self-similar patterns (Pietronero 1987; Ribeiro and Miguelote 
1998; Gabrielli et al. 2005), we have drawn straight lines for reference, whose slopes would 
approximately give the respective fractal dimensions of the distribution. The plots them- 
selves do not prove that the CNOC2 galaxy distribution is fractal for z > 0.1, but does 
not entirely rule out this possibility. In our view, this possibility clearly deserves further 
investigation. 

Of course, this is not saying that the spatial distribution of galaxies may be fractal. 
That is fitted by an EdS, or other FLRW, model and therefore approximately homogeneous 
on large scales. What it is saying is that the observationally determined average density of 
galaxies down our past light cone is possibly fractal. 

3.3. FLRW Cosmology with fi mo = 0.3 and Vt Ko = 0.7 

The differential densities can only be obtained in this cosmological model by means of 
numerical integrations. Here we shall not provide details of these calculations, which are 
the subject of a forthcoming paper (Iribarrem, Ribeiro and Stoeger 2006, in preparation), 
but limit ourselves in presenting data and plots containing similar results as in the EdS 
case. Suffice it to say that the methodology is precisely the same as described in §3.2: take 
the data in table 4 and apply them to equations (24) and (25), but evaluating numerically 
the geometrical part involving do, S and V in this FLRW cosmology. The results for all 
differential densities and distances are shown in table 6. 

Figures 5 and 6 show both differential densities against the redshift using the theoretical 
differential number count dN/dz obtained with the model, whereas figures 7 and 8 replace 
dN/dz with the observed [dN/dz] of table 4 obtained from the CNOC2 survey. One can 
clearly notice that all densities are strongly dependent on the distance choice, not only the 
theoretical, but the observed ones. However, the densities calculated with the galaxy area 
distance d G change from a small decrease in the EdS model to a more pronounced one in 
figures 7 and 8. That decrease does not happen in the theoretical plots shown in figures 5 
and 6, meaning that it can only be attributed to the CNOC2 data. One can also notice a 
clear power law pattern for densities plotted with the luminosity and redshift distances d L 
and d z at ranges where the redshift is higher than 0.1. 

Finally, figures 9 and 10 show both differential densities respectively plotted against d L 
and d z . As in the EdS model, there is a power law pattern in the tail of the plot, that is, for 
z > 0.1. That pattern seems to be even more linear than for the EdS cosmology. Straight 
lines for reference have also been drawn, showing the approximate slope of the linear decaying 
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behavior. Since power law patterns are usually related to self-similar fractal structures, the 
possibility that some standard cosmological models possess observational fractal patterns at 
ranges of z > 0.1 deserves further investigation in other galaxy survey samples. 

4. Conclusions 

In this paper we have extended the theory developed by Ribeiro and Stoeger (2003) for 
linking the relativistic cosmology source count to the luminosity function (LF) in order to 
extract the differential number count alN/dz from LF data. Since LF parameters assume a 
cosmological model, we developed the analytical means to factor out the assumed cosmology 
from the LF parameters, while preserving the LF analytical assumptions and source evolution 
corrections. We have applied this methodology to the CNOC2 galaxy redshift survey (Lin 
et al. 1999) and successfully extracted the observed differential number count [dN/dz] from 
this galaxy catalogue. Then we used this extracted data to carry out a statistical analysis 
based on the radial differential densities discussed in Ribeiro (2005). We have done so in 
the Einstein-de Sitter cosmology and the currently favored f2 mo = 0.3, £7a = 0.7 standard 
cosmological model. Our main aim was to test Ribeiro's (2005) proposition that different 
distance definitions when applied to galaxy survey analyses introduce ambiguities that render 
the statistical results unavoidably dependent on the distance choice. 

Our results show that Ribeiro's (2005) proposition is sound for both models under study. 
We calculated two types of radial statistical tools, namely the differential density 7 and the 
integral differential density 7*, with the CNOC2 derived data and found out that both have 
a strong dependency on the distance choice, as well as, showing an observational power law 
pattern for z > 0.1 when 7 and 7* are calculated with the luminosity distance d L and the 
redshift distance d z . 10 This power law decay is present in both cosmological models, EdS 
and fl mo = 0.3, Q\ = 0.7 standard model, although the linear power law tendency is more 
pronounced in the latter cosmology. If we cautiously interpret this power law pattern as 
possible evidence of a self-similar fractal structure, that would mean observational galaxy 
distribution fractal dimensions in the range D»1.6- 1.8 at scales of z > 0.1. 

In addition to the remarks of §3.1.2 above, it is helpful to remember that the cumu- 
lative number counts N is a dimensionless observational quantity with no dependence on 



10 Notice that this power law pattern only appears with some observational quantities. As stressed in 
§3.1.2 above, the results of this paper only refer to observational homogeneity as the underlying models 
still keep their spatial homogeneity, which is there by construction. See also Ribeiro (2001) for in depth 
discussions of this topic. 
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any volume definition. It can be used in other contexts, like star counting, in a purely 
observational setting, having no relation to theory. However, in determining related param- 
eters - (no) for example - we must know over what volume N has been counted, and this 
involves determining or assuming a definite cosmological model. In a cosmological context 
the theoretical counterpart of Nq, that is N, scales with the redshift according to the over- 
all mass-energy distribution of the universe, and this theoretical dependence given by the 
cosmological models is what we seek to relate to number count observations, that is, to N . 

The link between N and Nq, or between derivative differential number counts dN/dz 
and [dN/dz} , is provided by equation (1). In this equation dN is written in terms of the 
proper number density n. Since in equation (1) n is written in terms of proper volumes, it 
was necessary to use the area distance - to tie together the locally observed cross-sectional 
area, which enters in the proper volume, with the observed solid angle. So, we can only use 
equation (1) to obtain dN/dz from a number density if we write n in terms of the proper 
density. If n is not given in terms of the proper volume, a conversion is required, that being 
the case when one deals with the LF in its usual form. In theoretical calculations we can 
take the local mass-energy density p that appears on the right hand side of Einstein's field 
equation and get n from it (RS03, eq. 16). Then N = N(z) is independent of any volume, but 
still dependent on the overall matter distribution in the universe, which, through General 
Relativity, is given by the spacetime geometry. In particular, the affine parameter y is 
required in deriving N and the dependence of y on z is given by the cosmology specifically 
adopted (e.g., eq. 12 above). 

As a result of the procedure outlined above, the final dN/ dz is dimensionless and volume 
independent, meaning that we are then free to replace dN/dz as given by equation (1) with 
our observed, and also dimensionless, [dN/dz] Q in order to define other number densities in 
different volumes. The 7 and 7q functions are just different number densities, built in such a 
way as to show not only the scaling with z (dependence with the overall matter distribution), 
but also its dependence on one of the possible volume definitions of the assumed cosmological 
model (geometrical dependence). The idea is to make explicit the geometrical dependence 
on distance choice, whereas this choice has been implicit in most discussions about the 
large-scale galaxy distribution in the Universe. 

We should also mention that linking the observational counts and densities with spatial 
ones, and with spatial homogeneity, takes some work. One way one can do this is by 
using [dN/dz] data and Ellis' formula (1), or its observational version (eq. 4) (Stoeger 1987; 
Araujo and Stoeger 1999) to obtain n(z) at different redshifts, by observationally determining 
d A = d A (z) for at least some galaxies in the sample, and see if this n as function of z agrees 
with the FLRW theoretical results. This would really give us the proper densities on space- 
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like slices, and enable us to see spatial homogeneity. This is quite different from integrating 
over large volumes down our past light cone to get observationally averaged densities, which 
convolve proper densities from a wide range of redshifts. 

As a final remark, since spatial and observational homogeneities are two different fea- 
tures of a cosmological model, as extensively discussed in §3.1.2 above, it should be clear 
by now that further work needs to be done in this area, including a clear treatment of 
how one observationally tests for spatial homogeneity using observed galaxy distributions 
and other data, without first demonstrating or assuming that the cosmological region under 
consideration is "almost FLRW." 
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A. Error Analysis 

In RS03 the selection function for Schechter's type LF was written in terms of the 
absolute magnitude M w in the filter W as follows (RS03, eq. 96), 

^ w (z) = (O.41nlO)0JO a4(1+Q) ( M *- M "O e xp _i °- 4 ( M - M ^) dM w , (Al) 

J —oo 

where the three parameters that characterize the LF are the space density of galaxies 0*, 
the asymptotic slope of the faint end a, and the luminosity scale given in terms of absolute 
magnitude M*. The CNOC2 survey incorporates a luminosity scale evolution such that M* 
was rewritten in terms of two other parameters M[ and Q (Lin et al. 1999, eq. 10; RS03, eq. 
63) leading equation (Al) to be rewritten as follows, 

ip w (z) = 0.4 In 10 0, (A2) 

x / w 0A(l+a)[M>-Q( Z - .3)-M w ] ^ (_ 1q0 a[m:-Q(z-0.3)-M w ] j & ^ 

J —oo 

Lin et al. (1999) provided uncertainties for these four parameters, Ml, a, 0*, Q in their tables 
1-3. These errors were propagated in equation (A2) by means of standard error propagation 
techniques in order to produce the uncertainties appearing in tables 1, 2, 3 of this paper. 
Numerical and algebraic computing scripts were used for this task. 

To obtain uncertainties for the observed differential number counts [dN/dz] shown in 
table 4, we have used equation (17) to further propagate the already calculated errors for the 
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selection functions. Similarly, by means of equation (22) it was straightforward to evaluate 
uncertainties for the differential density 7. These results are shown in tables 5 and 6. 

To evaluate uncertainties for the integral differential density 7* we proceeded as follows. 
From equation (26) it is clear that £7* = 5N/V. Since we can write N = J (dN/dz)dz, it is 
easy to find SN, and, therefore, uncertainties of the integral differential density yield, 

Equation (17) and the already tabulated errors of the observed differential number counts 
5[dN/dz] allowed us to obtain the values for Sjq. The results are shown in tables 5 and 6. 
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Table 1. CN0C2 selection function a results from RS03 plus errors evaluated here (see 

appendix A) in the R c band vs. redshift. 

0.05 0.0181 ± 0.0038 0.0120 ± 0.0040 0.0220 ± 0.0130 

0.12 0.0182 ± 0.0038 0.0121 ± 0.0041 0.0221 ± 0.0130 

0.25 0.0183 ± 0.0039 0.0123 ± 0.0041 0.0223 ± 0.0130 

0.4 0.0185 ± 0.0039 0.0125 ± 0.0042 0.0225 ± 0.0132 

0.55 0.0186 ± 0.0040 0.0127 ± 0.0044 0.0227 ± 0.0135 

0.75 0.0187 ± 0.0040 0.0130 ± 0.0045 0.0230 ± 0.0141 

0.9 0.0187 ± 0.0041 0.0132 ± 0.0046 0.0232 ± 0.0148 

1.0 0.0187 ± 0.0041 0.0133 ± 0.0047 0.0234 ± 0.0153 



a Subscript numbers denote the spectral type morphology 
adopted for the CNOC2 survey (see §2.4 above). Units are 
Mpc~ 3 . 
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Table 2. CN0C2 selection function results from RS03 plus errors evaluated here (see 

appendix A) in the U band vs. redshift. 

w r 2 w 

0.05 0.0180 ± 0.0031 0.0119 ± 0.0038 0.0284 ± 0.0134 

0.12 0.0182 ± 0.0031 0.0121 ± 0.0038 0.0289 ± 0.0135 

0.25 0.0186 ± 0.0032 0.0123 ± 0.0039 0.0297 ± 0.0139 

0.4 0.0189 ± 0.0033 0.0126 ± 0.0040 0.0307 ± 0.0145 

0.55 0.0192 ± 0.0033 0.0128 ± 0.0042 0.0317 ± 0.0152 

0.75 0.0194 ± 0.0034 0.0131 ± 0.0043 0.0330 ± 0.0162 

0.9 0.0195 ± 0.0034 0.0134 ± 0.0045 0.0341 ± 0.0172 

1.0 0.0196 ± 0.0034 0.0135 ± 0.0045 0.0348 ± 0.0179 
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Table 3. CN0C2 selection function results from RS03 plus errors evaluated here (see 

appendix A) in the B AB band vs. redshift. 



B 



Z 4>2 $ 

0.05 0.0182 ± 0.0033 0.0122 ± 0.0036 0.0252 ± 0.0130 

0.12 0.0183 ± 0.0033 0.0123 ± 0.0036 0.0254 ± 0.0130 

0.25 0.0185 ± 0.0034 0.0125 ± 0.0037 0.0257 ± 0.0131 

0.4 0.0187 ± 0.0034 0.0128 ± 0.0038 0.0260 ± 0.0133 

0.55 0.0189 ± 0.0035 0.0130 ± 0.0040 0.0264 ± 0.0136 

0.75 0.0190 ± 0.0035 0.0133 ± 0.0041 0.0268 ± 0.0143 

0.9 0.0191 ± 0.0035 0.0135 ± 0.0042 0.0272 ± 0.0150 

1.0 0.0191 ± 0.0035 0.0136 ± 0.0043 0.0274 ± 0.0155 
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Table 4. CN0C2 differential counts for all wave bands and morphologies vs. redshift. 



dN 

dz 



0.05 


(2.914 ± 1.024) 


xlO 7 


0.12 


(1.394 ± 0.491) 


xlO 8 


0.25 


(4.409 ± 1.564) 


xlO 8 


0.4 


(8.125 ± 2.920) 


xlO 8 


0.55 


(1.142 ± 0.417) 


xlO 9 


0.75 


(1.487 ± 0.555) 


xlO 9 


0.9 


(1.679 ± 0.639) 


xlO 9 


1.0 


(1.780 ± 0.686) 


xlO 9 



Table 5. Redshifts, distances and differential densities calculated with CN0C2 differential 
number count data in an Einstein-de Sitter cosmology. Distances are in Mpc, densities in 

Mpc -3 and H = 100 km s _1 Mpc -1 . The differential densities were calculated using 
equations (24) and (25), where the observed differential number count [dN/dz] was given 
by the survey data results shown in table 4, whereas the second factor in the right hand side 
of equation (24) was obtained theoretically, that is, from the assumed cosmological model. 
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Table 6. Redshifts, distances and differential densities calculated with CN0C2 differential 
number count data in a Q mo = 0.3, Qa — 0.7 standard cosmology. Distances are in Mpc, 
densities in Mpc~ 3 and Hq = 100 km s _1 Mpc -1 . The differential densities were calculated 

similarly as explained in the caption of table 5. 
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Fig. 1. — Differential densities 70 against redshift in the EdS model with CNOC2 data. The 
same theoretical behavior of the differential count presented in Ribeiro (2005) is preserved 
with the CNOC2 differential number counting data. Notice the power law pattern when 
z > 0.1 for densities obtained with d L and d z . 
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Fig. 2. — Integral differential densities 7q against redshift in the EdS model. As in the 
previous plot, the critical dependence on the choice of distance is preserved in the CNOC2 
data, as well as the visible power law decaying pattern occurring when one chooses d L and 
d z as distances. 
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Fig. 3. — Plot of the differential densities obtained from the CNOC2 survey against the 
luminosity distance in EdS cosmology. Both densities decrease at higher d L and show a clear 
power law behavior in the tail. The line is for reference only, but its slope is such that in 
the language of fractals it would mean a fractal dimension D = 1.7. This plot does not 
prove that the galaxy distribution as presented in the CNOC2 survey data follows a fractal 
pattern, but that such an observationally self-similar structure is possible at least in some 
scales. 
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Fig. 4. — Plot of the differential densities obtained from the CNOC2 survey against the 
redshift distance in EdS cosmology. Similarly to the previous graph, both densities tend to 
follow a power law decay at higher distances. The straight line is for reference only, but 
its slope is such that it would mean D = 1.8. Again, one must remain cautious regarding 
interpreting this self-similar behavior in the tail as evidence of a fractal pattern. However, 
the plot indicates that this possibility should not be ruled out either. Further investigation 
with other samples is clearly required to clarify this point. 
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Fig. 5. — Plot of the theoretical differential density 7 for each observational distance in the 
fi mo = 0.3, J1a — 0-7 standard cosmology. The results are similar to those of EdS model 
as shown in Ribeiro (2005) and are presented here for comparison with the observationally 
derived 70 shown in figure 7. Notice that similarly to the EdS case 7 A blows up at about 
z ~ 1.5, but that in itself does not mean that the volume defined with d A is incorrect, but 
only that the way 7 A is defined has a singularity, which has no special physical significance 
other than that the distance it involves has a critical point. 
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Fig. 6. — Plot of the theoretical integral differential density 7* for each observational distance 
in the f2 mo = 0.3, Q\ = 0.7 standard cosmological model. The results are similar to those 
of EdS cosmology (Ribeiro 2005). They are presented here so that one can compare with 
the observationally derived 7o plotted in figure 8. 
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Fig. 7. — This figure shows the observed differential density 70 vs. the redshift in the 
H mo = 0.3, Qa = 0.7 standard cosmology with CNOC2 differential number counts data of 
table 4. Notice that similarly to the results shown in figure 1 there is as power law pat- 
tern at the tail of the plot for both d L and d z . However, the differential density calculated 
with the galaxy area distance d G does not remain unchanged, but also decays at higher z, 
an effect that can only be explained because this plot uses the CNOC2 differential number 
count [dN/dz] for calculating 70, rather than its theoretically derived values used in figure 5. 
Therefore, it is possible that either the CNOC2 survey is undercounting galaxies at larger z 
or that its LF source evolution equation is not correct. A third, less likely, possibility is that 
there could be real spatial inhomogeneity and in this case the FLRW model would not be 
correct. 
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Fig. 8. — Plot of the observed integral differential density 7q against the redshift in the 
H mo = 0.3, Q\ = 0.7 standard cosmology with a CNOC2 LF derived [dN/dz] Q . As in the 
previous plot one can also notice a power law pattern for both the luminosity distance d L 
and redshift distance d z . But, still as in the previous plot, 7q calculated with d G decreases 
at higher z, departing from the theoretical behavior presented in figure 6. 
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Fig. 9. — Differential densities against the luminosity distance in the Q mo = 0.3, ^a = 0.7 
standard cosmology using the CNOC2 differential number count data of table 4. The power 
law behavior is visible in the z > 0.1 ranges. The straight line is for reference only but its 
slope would mean a fractal dimension of approximately D — 1.7. As cautioned above, one 
must remain careful in interpreting this pattern as evidence of a fractal structure, since the 
results were obtained only with the CNOC2 survey data. 
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Fig. 10. — Differential densities versus the redshift in the f2 mo = 0.3, £l^ = 0.7 standard cos- 
mological model where [dN/dz} was derived from the LF obtained with the CNOC2 galaxy 
survey (see table 4). Similarly to the previous figure, one can notice a power law pattern in 
the tail of both plots. The straight line is for reference only, but would mean approximately 
D — 1.8 if one interprets these results as possible evidence of a fractal structure for z > 0.1 
ranges. 



